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Abstract. We produce a 10-day series of simulated Dopplcr images at a 
15-minute cadence that reproduces the spatial and temporal characteristics seen 
in the SOHO/MDI Doppler data. Our simulated data contains a spectrum of 
cellular flows with but two necessary components — a granule component that 
peaks at wavenumbers of about 4000 and a supergranule component that peaks 
at wavenumbers of about 110. We include the advection of these cellular compo- 
nents by a differential rotation profile that depends on latitude and wavenumber 
(depth). We further mimic the evolution of the cellular pattern by introducing 
random variations to the amplitudes and phases of the spectral components at 
rates that reproduce the level of cross-correlation as a function of time and lat- 
itude. Our simulated data do not include any wave-like characteristics for the 
supergranuies yet can accurately reproduce the rotation characteristics previ- 
ously attributed to wave-like characteristics. 



1. Introduction 

Supergranuies are cellular flow structures observed in the solar photosphere with 
typical diameters of about 30 Mm and lifetimes of about one day. They cover 
the entire surface of the Sun and are intimately involved with the structure and 
evolution of the magnetic field in the photosphere. The magnetic structures of 
the chromospheric network form at the boundaries of these cells and magnetic 
elements are shuffled about the surface as the cells evolve. The diffusion of 
magnetic elements by the evolving super granuies has long been associated with 
the evolution of the Sun's mag netic fie ld ILeightonl ([1964ft ] . 



Supergranuies were discovered by iHartJ ( 19541b While these cellular flow s 



were quickly identified as convective features Leighton. Noves. Sz Simon (1962)] 
the difficulty of detecting any associated thermal features consistent with that 
dentifica tion (i.e . hot c ell centers) has made this identification somewhat prob- 
lematic [Wordenl ([1975ft ]. 



T he rota t ion of the supergranuies has added further mystery to their na- 
ture. iDuvalll (H980) cross-correlated the equatorial Doppler velocity patterns 



and found that the supergranuies rotated more rapidly than the plasma at the 
photosphere and that even faster rotation rates were obtained when longer (24- 
hour vs. 8-hour) time interval s were used. He attrib uted this be h avior to a 
surface shear layer [proposed bv lFoukal fc JokipiH ( 19751 ) and iFoukal ( 19791 ) and 
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modeled by iGilman &: Foukall ( 19791 )] in which larger, longer-lived, cells e xtend 
deeper into the more rapidly rotating layers. ISnodgrass & Ulrichl (|l99d ) used 
data from Mount Wilson Observatory to find the rotation rate at different lati- 
tudes and noted that the rotation rates for the Doppler pattern were some 4% 
faster than the spectroscopic rate and, more interes tingly, some 2% faste r than 
the magnetic features and sunspots. More recently iBeck &: Schoul (|2000i ) used 
a 2D Fourier transform method to find that the Doppler pattern rotates more 
rapidly than the shear layer itself and that larger features do rotate more rapidly 
than the smaller features. They suggested that supergranules have wave-like 
characteristics with a pr eference for prograde propagati on. 

In a previous paper [Hathawav. Williams. & Cuntd (|2006T )] we showed that 
this "super-rotation" of the Doppler pattern could be attributed to projection 
effects associated with the Doppler signal itself. As the velocity pattern ro- 
tates across the field of view its line-of-sight component is modulated in a way 
that essentially adds another half wave and gives a higher rotation rate that 
is a function of wavenumber. In that paper we took a fixed velocity pattern 
(which had spatial characteristics that matched the SOHO/MDI data) and ro- 
tated it rigidly to show this "super-rotation" effect. While this indicated that 
this Doppler projection effect should be accou nted for, the fixed pa ttern could 
not account for al l the v ariations reported by IBeck Sz Schoul ( 200dl ) . Further- 
more, when lSchou ( 20031 ) "divided-out" the line-of-sight modulation he still saw 
prograde and retrograde moving components. 

In this paper we report on our analyses of simulated data in which the su- 
pergranules are advected by a differential rotation that varies with both latitude 
and depth. The data is designe d to fa ithful ly mim i c the SOHO/MDI data that 
was analyzed in IBeck &: Schoul ( 20001 ) and ISchoul ( 20031 ) and the analyses are 
reproductions of those done in earlier studies. 



The Data 



The full-disk Doppler images from SOHO/MDI [Scherrer et all (|l995l )] are ob 



tained at a 1-minute cadence to res olve the tempora l variat ions associated with 
the p -mode oscillations. We [c.f. lHathawav et al.l (200(3) and IBeck Sz Schou 



(2000)] have temporally filtered the images to remove the p-mode signal by us- 



ing a 31-minute long tapered Gaussian with a FWHM of 16 minutes on sets 
of 31 images that were de-rotated to register each to the central image. Series 
of these filtered images were formed at a 15-minute cadence over the 60-day 
MDI Dynamics Runs in 1996 and 1997. This filtering process effectively re- 
moves the p-mode signal and leaves behind the Doppler signal from flows with 
temporal variations longer than about 16 minutes. Supergranules, with typical 
wavenumbers of about 110, are very well resolved in this data (at disk center 
wavenumbers of about 1500 are resolved). While granules are not well resolved, 
they do appear in the data as pixel-to-pixel and image-to-image "noise," as a 
convective blue shift (due to the correlation between brightness and updrafts), 

and as resolved structures for the largest members. 

The simulated data are constructed in the manner described in Hathawav 



(|1988T ). iHathawavi (| 19921 ). lHathawav etaD ((2000), and lHathawav et al.1 (12002) 
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from vector velocities generated by an input spectrum of complex spectral coef- 
ficients for the radial, poloidal, and toroidal components. 

To simulate the observed line-of-sight velocity the three vector velocity com- 
ponents are calculated on a grid with 1024 points in latitude and 2048 points 
in longitude. A Doppler velocity image is constructed by determining the lon- 
gitude and latitude at a point on the image, finding the vector velocity at that 
point using bi-cubic interpolation, and then projecting that vector velocity onto 
the line-of-sight. The line-of-sight velocities at an array of 16 points within each 
pixel are determined and the average taken to simulate the integration over a 
pixel in the acquisition of the actual MDI Doppler data. 

With the current simulations we have added two changes that were not in- 
cluded in our previous work on individual Doppler images. First, we add velocity 
"noise" at each pixel. This represents the contribution from the spatially un- 
resolved granules that, nonetheless, have temporal varibility that is not filtered 
out by the 31-minute temporal filter. This noise has a center-to-limb variation 
due to the foreshortened area covered by each pixel and is randomly varied from 
pixel to pixel and from one Doppler image to the next. The noise level is de- 
termined by matching the initial drop in correlation from one image to the next 
that is seen in the MDI data. Secondly, we treat the instrumental blurring in 
a more realistic manner. Previously we took the Doppler velocity image and 
convolved it with an MDI point-spread-function. We now make red and blue 
intensity images from our Doppler velocity image and a simple limb darkened 
intensity image, convolve those with an MDI point-spread-function, and con- 
struct a blurred Doppler velocity image from the difference divided by the sum. 
This process yields a Doppler velocity image that is virtually indistinguishable 
from an MDI Doppler velocity image. 

The velocity pattern is evolved in time by introducing changes to the spec- 
tral coefficients based on two processes - the advection by an axisymmetric zonal 
flow (differential rotation) and random processes that lead to the finite lifetime 
of the cells. 

The advection by the differential rotation is governed by an advection equa- 
tion 



where u is a vector velocity component and is the differential rotation 

profile. We represent u as a series of spherical harmonic components and project 
this advection equation onto a single spherical harmonic which gives a series 
of coupled equations for the evolution of the spectral coefficients. Solid body 
rotation simply introduces a phase variation for each coefficient. Differential 
rotation couples the change in one spectral coefficient to spectral coefficients 
with wavenumbers £ ± 2 and £ ± 4 for differential rotation dependent on cos 2 9 
and cos 4 8. 

The finite lifetimes for the cells are simulated by introducing random per- 
turbations to the spectral coefficient amplitudes and phases. The size of these 
perturbations increases with wavenumber to give shorter lifetimes to smaller 
cells. 



du 
~dt 




(1) 
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3. The Analyses 



Several anaylsis programs were applied to both the MDI data and the simulated 
data. Convectio n spectra for indi vidua l images w e re ob tained using the meth- 
ods described by lHathawavl ( 19871 ) and lHathawavl ( 19921 ) — the Doppler signal 
due to the motion of the observer is removed, the convective blue shift signal 
is identified and removed, the data is mapped to heliographic coordinates, the 
axisymmetric flow signals due to differential rotation and meridional circulation 
are identified and removed, and the remaining signal is projected onto spherical 
harmonics. The averaged spectra from the 1996 MDI Dynamics Run and from 
our 10-day simulated data run are nearly perfectly matched at all wavenum- 
bers. This match is obtained by adjusting the input spectrum for the simulated 
data. This spectrum contains two Lorentzian-like spectral components — a su- 
pergranule component centered on £ ~ 110 with a width of about 100 and a 
granule component centered on £ ~ 4000 with a width of about 4000. The MDI 
spectrum is well matched wi th just these two comp onents without the addition 
of a mesogranule component November et al. ( 198ll )]. In fact, we find a distinct 
dip in the spectrum at wavenumbers £ ~ 500 that should be representative of 
mesogranules. This dip is also seen in spectra of the MDI high resolution data 



Hathaway et all (|2000)]. 



Additional analyses are applied to the data after it has been mapped onto 
heliographic coordinates. Longitudal strips of this data centered on latitudes 
from 75° south to 75° north w ere cros s -corre late d with corresponding strip s 
from later images as was done by iDuvalj ( 19801 ) and lSnodgrass fc Ulrichl (j 1990I K 
The longitudinal shift of the cross-correlation peak gives the rotation rate while 
the height of this peak is associated with cell lifetimes. These strips were also 
Fourier analyzed in longitude to get spectral coefficients and t hose coefficents 
were F ourier anaylzed in time over 10-day intervals as was done bv lBeck Sz Schoul 
(2000) to get rotation rates as functions of wavenumber. 

Fig. 1 shows the rotation rates from the cross-correlation analysis. The 
rotation profiles from the simulated data match those from t he MDI data. Both 
sh ow faster rotation ra t es for longer time lags as noted by iDuvalll ([l980) and 
by ISnodgrass fe Ulrichl ( 19901 ). This indicates that we have found the right 
latitudinal differential rotation profile. 

The strength of the correlations as functions of latitude and time lag for 
both the MDI data and the simulated data are also well matched. This indicates 
that we have found the right lifetimes for the c ells. 

We have also reproduced the analysis of ISchoul ( 20031 ). The data strips 
are apodized and multiplied by longitude dependent functions designed to re- 
move the Doppler projection effect and to isolate either longitudinal motions or 
latitudinal motions. The strips are shifted in longitude according to the differ- 
ential rotation rate and then Fourier analysed in space and time to obtain "kcu" 
diagrams. 

Fig. 2A shows the equatorial rotation rate as a function of wavenumber 
for the simulated data wh ile Fig. 2B sh o ws th e ku diagram. These should 
be compared to Fig. 4 of iBeck fc School (|200d > and Fig. 8 of ISchoul (|2003l ) 
respectively. 
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Figure 1. Rotation profiles from cross-correlation analyses of MDI data 
(solid lines) and simulated data (dotted lines) . All profiles match at virtually 
all latitudes and time lags (Panel A — 2-hour, Panel B — 4-hour, Panel C 
- 8- hour, Panel D — 1 6- hour). Note that the measured equatorial rotation 
rate increases with time lag for both datasets. 




Figure 2. Panel A. The equatorial rotation rate as a function of wavenum- 
ber from our analysis of the simulated data showing the increase at small 
wavenumbers. Most of this increase is due to the Doppler velocity projection. 
Panel B. The ku diagram from our analysis of the simulated data showing 
both prograde and retrograde moving elements. This spread in power is due 
to the random evolution of the cellular pattern. 
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4. Conclusions 

We have produced simulated data in which the cellular structures (supergran- 
ules) are advected by differential rotation and evolve by uncorrelated random 
changes. When we compare results from analyses of this data with those from 
analyses of the MDI data we find that the simulated data exhibits the same 
characteristics as the MDI data — the visual structures, the power spectra, the 
rotation characteristics, and the evolution rates all match. Whil e some of these 



charac teris tics have bee n attributed to wave-like properties [c.f. iBeck Sz Schou 



(200(3) andlSchoul (120031 )] our simulated data is simply advected by a zonal flow 



(differential ro tation) with s peeds that never exceed those determined from he- 
lioseismology [Schou et al.l (|l998l )]. The differential rotation we impose does, 
however, have a dependence on wavenumber I. If we assume that the rotation 
rate of cells with diameters, D = 2ttRq/£, reflects the rotation rate at a depth, 
d = D/2, then the surface shear layer indicated by our differential rotation has 
a thickness of about 20 Mm — somewhat thinner than the 30 Mm suggested by 
helioseismic inversions [Schou et al.l (|l998l )]. 
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